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We report a ^Li nuclear magnetic resonance (NMR) study of a frustrated three-dimensional spin-| 
antiferromagnet Li2CuW208 and also explore the effect of non-magnetic dilution. The magnetic 
long-range ordering in the parent compound at Tn — 3.9 K was detected from the drastic line 
broadening and a peak in the spin-lattice relaxation rate (1/Ti). The NMR spectrum above Tn 
broadens systematically, and its full width at half maximum (FWHM) tracks the static spin sus¬ 
ceptibility. From the analysis of FWHM vs. static susceptibility, the coupling between the Li nuclei 
and Cu^^ ions was found to be purely dipolar in nature. The magnitude of the maximum exchange 
coupling constant is Tmax/^B — 13 K. NMR spectra below Tn broaden abruptly and transform into 
a double-horn pattern reflecting the commensurate nature of the spin structure in the ordered state. 
Below Tn, l/Ti follows a T® behavior. The frustrated nature of the compound is confirmed by per¬ 
sistent magnetic correlations at high temperatures well above Tn. The dilution of the spin lattice 
with non-magnetic Zn atoms has dramatic influence on Tn that decreases exponentially similar to 
quasi-one-dimensional antiferromagnets, even though Li 2 CuW 2 08 has only a weak one-dimensional 
anisotropy. Heat capacity of doped samples follows power law (Cp oc r“) below Tn, and the expo¬ 
nent (a) decreases from 3 in the parent compound to 1 in the 25% doped sample. 

PACS numbers: 75.10.Jm, 75.30.Et, 75.50.Ee, 76.60.-k, 75.30.Kz 


I. INTRODUCTION 

Geometrically frustrated spin-1/2 triangular systems 
have attracted considerable attention due to numerous 
quantum phenomena that they exhibit at low temper¬ 
atures. The combination of geometric frustration, low 
spin value, and reduced dimensionality significantly en¬ 
hances quantum fluctuations which can lead to the for¬ 
mation of exotic ground states.^ The most striking ex¬ 
ample is the formation of a quantum spin liquid (QSL), 
which is a strongly correlated disordered ground state. 
Anderson suggested that the spin-1/2 Heisenberg trian¬ 
gular antiferromagnets have a resonating-valance-bond 
(RVB) ground state, one kind of the QSL.^ Indeed, sev¬ 
eral triangular antiferromagnets, such as the spin-1 com¬ 
pound NiGa 2 S 4 and spin-1/2 organic materials like k- 
(BEDT-TTF) 2 Cu 2 (CN )3 and EtMe3Sb[Pd(dmit)2]2, do 
not show long-range magnetic order (LRO) and are even 
established as QSL candidates.On the other hand, the 
majority of triangular antiferromagnets entail subtle de¬ 
viations from the ideal regime that may reduce the frus¬ 
tration and stabilize the LRO. In fact, even the idealized 
case of the quantum Heisenberg model on the triangular 
lattice features LRO with the 120° spin structure®’^ and 
thus lacks the spin-liquid phase. Most of the spin-i tri¬ 
angular antiferromagnets are magnetically ordered at low 
temperatures. They may also show plateaux at 1 /3 of the 
saturation magnetization, where the system evolves from 
a 120° magnetic ordering phase to an up — up — down 
{uud) phase in a finite field,® as in spin-i compounds 


Cs 2 CuBr 4 and Ba 3 CoSb 209 .® 

Li 2 CuW 208 crystallizes in a triclinic structure with 
space group PI (Ref. 12) and features a unique spin lat¬ 
tice, which is strongly frustrated along all three crystal¬ 
lographic directions. The interaction topology in the ab 
plane resembles a triangular lattice, whereas interplane 
couplings are oblique and likewise frustrated.^® There is 
one Cu site, two equivalent W sites, and two equivalent Li 
sites in the unit cell. The CuOg octahedra that can also 
be viewed as Cu 04 planar units are connected through 
WOe octahedra. The Li atoms are located close to each 
of the Cu^“'" triangular layer. Figure 1 shows the trian¬ 
gular layer formed by Cu^“'" ions and the collinear ar¬ 
rangement of spins inferred from the neutron diffraction 
experiments.^® The Li atom is coupled to two neighbour¬ 
ing Cu^’*' spins (either up — up or down — down) of the 
same plane. 

Recently, we reported magnetic susceptibility x(T), 
heat capacity Cp{T), and neutron diffraction on 
Li2CuW208.®® Clear signatures of the magnetic three- 
dimensionality and strong frustration were inferred from 
these measurements and rationalized microscopically us¬ 
ing band-structure calculations, x(^) s-ud Cp{T) exhibit 
broad maxima at around 8.5 K and 6 K, respectively, 
associated with the short-range order and then a kink 
(x) or a A-type anomaly (Cp) at Tn ~ 3.9 K associated 
with the magnetic long-range ordering (LRO). Neutron 
diffraction experiments reveal collinear magnetic order 
below 3.9 K with the propagation vector (0,^,0) and the 
ordered moment of 0.65(4) /xb, which is reminiscent of a 
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FIG. 1. The triangular layer formed by ions, columnar 

arrangement of spins, and the couplings of Li atoms to the 
magnetic spins are shown. 


2D antiferromagnet even though Li2CuW208 is clearly 
three-dimensional (3D). Remarkably, the collinear order 
is not anticipated on the classical level, where an in¬ 
commensurate spiral structure would have lower energy. 
Therefore, a strong stabilization of the collinear order by 
quantum fluctuations has been envisaged. 

Li2CuW208 is isostructural to its Co and Ni ana¬ 
logues. Both Li2NiW208 and Li2CoW208 compounds 
undergo two successive magnetic transitions at low tem¬ 
peratures without having any short-range order above 
Many of the triangular-lattice compounds ex¬ 
hibit two phase transitions at low temperatures too.^® 
The presence of only one magnetic transition and the ro¬ 
bust collinear order in Li 2 CuW 208 are clearly unusual 
and require detailed investigation. 

In this paper, we report ^Li NMR measurements on 
Li 2 CuW 2 08 and the effect of non-magnetic Zn^+ doping 
at the magnetic Cu^+ site. Compared to the previous 
study where only bulk measurements were performed, 
we use NMR, which is a powerful local tool for studying 
magnetic interactions and magnetic order in frustrated 
spin systems. The reduction in the ordered magnetic mo¬ 
ment (0.65 ^b) below its classical value of I has been 
ascribed to strong quantum fluctuations, but alterna¬ 
tively it could be understood as a coexistence of ordered 
and disordered phases below Tjv. Our NMR study rules 
out this possibility and corroborates the formation of the 
commensurate order throughout the sample. We also 
demonstrate the persistence of strong spin correlations 
well above T^r, which is a hallmark of strongly frustrated 
magnets, and analyze spin dynamics of Li 2 CuW 208 . 
The experiments on magnetic dilution via Zn^"*' dop¬ 
ing bring yet another perspective on the frustrated mag¬ 
netism of Li2CuW208. We demonstrate an abrupt de¬ 
crease in the Neel temperature upon doping, an effect 
reminiscent of quasi-one-dimensional (ID) antiferromag- 


nets and very peculiar considering the largely 3D nature 
of Li 2 CuW 2 08 . 

II. EXPERIMENTAL 

Poly crystalline samples of Li 2 (Cui_ 2 :Zna;)W 208 (x = 
0, 0.05, 0.10, 0.15, 0.20, 0.25, and 0.30) were prepared 
by the conventional solid-state reaction technique using 
Li 2 C 03 (Aldich, 99.999%), CuO (Aldich, 99.999%), ZnO 
(Aldich, 99.99%), and WO 3 (Aldich, 99.999%) as start¬ 
ing materials. The stoichiometric mixtures were ground 
and fired at 650° C for 24 hours to allow decarbonation 
and then at 700°C for 48 hours with one intermediate 
grinding. Phase purity of the samples were checked using 
powder x-ray diffraction (XRD) experiment (PANalyti- 
cal. Empyrean) with Cu Ka (A = 1.54060 A) radiation. 
The samples with x upto 25%, were single-phase, but 
at higher doping concentrations several impurity phases 
including Li 2 W 04 emerged. Our repeated attempts to 
achieve higher doping levels by increasing or lowering the 
firing temperature were unsuccessful. Therefore, we fo¬ 
cus on studying the samples with x < 30%, where a minor 
amount of non-magnetic impurities does not hinder the 
data analysis. 

Temperature (T) dependent magnetic susceptibility 
x(T) and heat capacity C'p(T) measurements were per¬ 
formed using a commercial Physical Property Mea¬ 
surement System (PPMS, Quantum Design). For the 
x(T) measurement, the vibrating sample magnetometer 
(VSM) attachment to the PPMS was used. C'p(T) was 
measured by the relaxation technique on a pressed pellet 
using the heat capacity option of the PPMS. For the low 
temperature (T < 2 K) Cp{T) measurements, a Heliun-3 
attachment was used. 

Nuclear magnetic resonance (NMR) experiments were 
carried out using pulsed NMR technique on the ^Li nu¬ 
cleus with the nuclear spin 7 = 3/2 and gyromagnetic 
ratio xn/Stt = 16.546 MHz/T. The NMR spectra were 
obtained either by Fourier transform of the NMR echo 
signal at a fixed field of 77 = 1.5109 T or by sweep¬ 
ing the magnetic field at a corresponding fixed frequency 
of 24.79 MHz. ^Li spin-lattice relaxation rate l/Ti was 
measured using a conventional saturation pulse sequence. 
^Li spin-spin relaxation rate I/T 2 was obtained by mea¬ 
suring the decay of the echo integral with variable spacing 
between the tt/2 and tt pulses. 

Electronic structure and electric field gradient on the 
Li site were calculated in the FPLD code^^ using the same 
procedure as in Ref. 13. 

III. RESULTS 
A. Crystallography 

Le Bail fit of the powder XRD data was performed 
using the FullProf software package based on the tri- 
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FIG. 2. X-ray powder diffraction pattern (open circles) at 
room temperature for Li2CuW208. The solid red line repre¬ 
sents the calculated pattern using Le Bail fit, with the vertical 
bars showing the expected Bragg peaks, and the lower solid 
blue line representing the difference between the observed and 
calculated intensities. 


clinic structure with the space group PI to determine 
the lattice parameters.^® All the data sets could be fit¬ 
ted using the structural data from Ref. 12 of the parent 
compound as the initial parameters. Figure 2 shows the 
room-temperature powder XRD pattern of Li 2 CuW 2 08 
along with the calculated pattern. The obtained best fit 
parameters are listed in Table I for all the Zn doped sam¬ 
ples. For the parent compound, the obtained lattice pa¬ 
rameters are in close agreement with the literature data 
[a = 4.9669(1) A, b = 5.4969(1) A, c = 5.8883(1) A, 
a = 70.72(1)°, /3 = 85.99(1)°, and 7 = 66 . 04 ( 1 )°]. 

In the doped samples, the lattice parameter a is de¬ 
creasing while b and c are increasing with increasing x 
(Table I). Similarly, the angles a and 7 are decreasing 
while /3 is increasing with increasing x. The over all unit 
cell volume {V) increases systematically, but the change 
in V is marginal, in agreement with the fact that the 
ionic radius of Zn^+ ( 0.6 A) is only slightly larger than 
that of Cu^+ (0.57 A). 

B. ^Li NMR spectra 

For the quadrupolar ^Li (/ = 3/2) nucleus, one expects 
two satellite lines along with the central line correspond¬ 
ing to three allowed transitions. As shown in Fig. 3, 
we observed only a narrow single spectral line without 
any satellites, which is likely due to the low quadrupo¬ 
lar frequency or distribution of intensity of the satellites 
over a broad frequency range. This type of single spec¬ 
tral line corresponding to the Zeeman central transition 
(-1-1/2 ^1/2) is commonly observed in ^Li NMR on 
low-dimensional oxides.Indeed, we calculated elec¬ 
tric field gradient on the Li site and obtained a rather 
low value of 14^ = 8.89 x 10^®V/m^ that yields the 
quadrupolar frequency of 40 kHz, which is comparable 
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FIG. 3. ^Li Fourier transform spectra measured at an applied 
field of i7 = 1.5109 T and at different temperatures. 


to the linewidth. 

No shift in the central peak position was observed over 
the whole measured temperature range, above Tn- The 
NMR shift K is directly proportional to the spin sus¬ 
ceptibility (xspin) and can be written as K = -^Xspin, 
where the proportionality constant A is the transferred 
hyperfine coupling between the ^Li nucleus and neigh¬ 
boring Cu^'*' spins, and Na is the Avogadro’s number. 
Thus our experimental observation of the temperature- 
independent line shift reflects a weak transferred hyper¬ 
fine coupling of the ^Li nuclei with the Cu^“'" spins. 

The spectra were found to broaden progressively upon 
lowering the temperature. Below about 75 K (see Fig. 3), 
the line shape becomes anisotropic. Figure 4 shows the 
full width at half maximum (FWHM) or Av of the ^Li 
NMR spectra as a function of temperature, above Tn- 
Av is increasing with decreasing temperature and then 
shows a broad maximum at around 10 K similar to the 
x(T) data.^® For a comparison, in the same Fig. 4, x(T) 
is also plotted in the right y-axis. 

In a polycrystalline sample, the NMR line broadening 
originates from various sources. The principal sources 
are defects or disorder at the magnetic sites which induce 
staggered magnetization, anisotropy in x(T), the dipolar 
interaction between the nuclei, and the dipolar hyperfine 
interaction of the nuclei with the neighboring magnetic 
ions. A small variation in bond length between the nuclei 
and the magnetic ions also contribute to the line broad¬ 
ening. In Li 2 CuW 208 , since the transferred hyperfine 
coupling between ^Li nuclei and the neighboring Cu^®" 
ions is negligibly small, the ^Li NMR line width should 
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TABLE I. The refined lattice constants (a, b, and c), angles [a, /?, and 7 ), and the unit cell volume V obtained from the Le 
Bail fit of the powder XRD data for Li 2 (Cui_i,Zna;)W 208 . 


X a (A) b (A) c (A) 

0 4.9627(2) 5.4925(2) 5.8836(2) 
0.05 4.9596(1) 5.5069(1) 5.8892(2) 
0.1 4.9558(1) 5.5152(2) 5.8907(2) 
0.15 4.9524(2) 5.5218(4) 5.8910(2) 
0.2 4.9487(3) 5.5298(4) 5.8918(1) 
0.25 4.9434(2) 5.5353(2) 5.8935(1) 
0.3 4.9408(1) 5.5446(1) 5.8941 (1) 


«(°) pn 7 (°) 

70.717(1) 85.99(1) 66.04(1) 
70.5812(9) 86.179(1) 65.917(1) 
70.4937(8) 86.3250(9) 65.824(1) 
70.4180(7) 86.439(8) 65.7422(9) 
70.3331(8) 86.5800(1) 65.6466(2) 
70.230(1) 86.723(2) 65.55(2) 
70.147 (2) 86.846(2) 65.483(1) 


V (A^) 
160.3729(1) 
160.8522(4) 
161.0104(2) 
161.0986(8) 
161.2340(8) 
161.2703(1) 
161.4730(1) 



1.5 

1.0 

0.5 

0.0 


FIG. 4. Temperature dependence of the full width at half 
maximum (Aiy) of the ^Li NMR spectra measured at an ap¬ 
plied field of 77 = 1.5109 T (left y-axis) and magnetic suscep¬ 
tibility (x) measured at 77 = 1 T (right y-axis). Inset shows 
the Ai/ vs. X with temperature as an implicit parameter and 
the solid line is the fit by Eq. (3). 


not be very sensitive to the staggered magnetization and 
the anisotropy in x{T). Similarly, if there are distor¬ 
tions at the lattice site, especially in the powder sample, 
this leads to a small variation of inter-atomic distances. 
From our temperature dependent x-ray diffraction and 
neutron scattering experiments, no significant change in 
bond length between Li and Cu^+ ions was found as a 
function of temperature. If at all, a small variation is 
there, it will not have any pronounced effect on the line 
width and hence ignored. 

Therefore, we presumed that the shape and width of 
our NMR spectra are a measure of two main interactions: 
(i) the dipolar interaction between the nuclei and (ii) the 
dipolar hyperfine interaction of the ^Li nuclei with the 
neighboring Cu^"*" magnetic ions. So, one can write Av 
above Tn as^^ 


interaction between the ^Li nuclei and magnetic Cu^"*" 
ions. In a polycrystalline sample, under magnetic field, 
this anisotropic coupling give rise to inhomogeneous line 
broadening. It is proportional to the bulk magnetic sus¬ 
ceptibility (x) and can be written as 


(^^)aniso 

V-L 




< ^ > 
r^H ’ 


( 2 ) 


where vl = denotes the ^Li Larmor frequency for an 
applied external magnetic field H, is the average dipo¬ 
lar coupling constant between the Li nucleus and mag¬ 
netic Cu^+ ions, and r is the average distance between 
them. 

It is to be noted that the shape of our experimental 
spectra doesn’t appear to be due to what one would ex¬ 
pect exactly from the anisotropy of the dipolar interac¬ 
tion. Perhaps, there are some other sources which also 
partly contribute to the line broadening. 

Putting Eq. (2) in Eq. (1) we can write 


Ah' = Avq + (3) 

In the inset of Fig. 4, Ai/ is plotted as a function of x with 
temperature as an implicit parameter. Av vs. x follows 
a linear behaviour and is fitted by Eq. (3) in the temper¬ 
ature range 9 K< T <175 K. The parameters obtained 
from the fit are Avq — 18.4 KHz and A^ ~ 7.83x10^^ 
cm“^. This value of A^ is similar to the values reported 
for several other compounds^^ and is of the right order 
of magnitude with the expected dipolar interaction of Li 
nuclei with Cu^’*' ions at a mean distance of 2.93 A apart. 
Thus we conclude that the hyperfine interaction of the Li 
nucleus with Cu^"*" moments is largely of dipolar origin. 


C. ^Li relaxation rates (l/7i and 1 / 72 ) 


Av = Avo + {Av) 

aniso; ( 1 ) 

where A^'o is the the contribution from the dipolar in¬ 
teraction between the nuclei, which is temperature- and 
field-independent. The second term (Aiz)aniso is the 
powder average over the anisotropic dipolar hyperfine 


The recovery of the longitudinal magnetization after a 
saturation pulse was fitted well by a single exponential 
function 


1 - 


M{t) 

Mo 




( 4 ) 
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FIG. 5. (a) Temperature dependence of the ^Li spin-lattice 
relaxation rate, l/Ti, for Li2CuW208 aX H — 1.5109 T. The 
solid line and the dotted line correspond to l/Ti oc T® and 
1/Ti oc T^g-A/feeT i30]2aviors, respectively. The dashed line 
is the linear fit (1/Ti oc T) in the 5.5 K < T < 16 K range. 
The inset shows l/Ti vs. —log// measured at T = 8 K. (b) 
l/yTiT as a function of temperature. 


where M(t) is the nuclear magnetization at a time t after 
the saturation pulse and Mq is the equilibrium value of 
the magnetization. The temperature dependence of 1/Ti 
extracted from the fitting is shown in Fig 5(a). At high 
temperatures (T > 20 K), 1/Ti is almost temperature- 
independent (in the paramagnetic regime), which is of¬ 
ten observed in a system with localized moments when 
the temperature is higher than the exchange energy be¬ 
tween the spins.With decreasing temperatures, 1/Ti 
decreases almost linearly for T < 20 K where the com¬ 
pound starts approaching the short-range-ordered regime 
(T ~ 10 K). With further decrease in temperature, a 
weak but discernible anomaly was observed at Ty — 
3.6 K due to the critical slowing down of fluctuating mo¬ 
ments while approaching the long-range magnetic order¬ 
ing, nearly consistent with the thermodynamic^^ mea¬ 
surements (Tn — 3.9 K).^^ Below Tn, 1/Ti decreases 
rapidly towards zero as a result of the disappearance of 
the critical fluctuations in the ordered state. 

We also measured 1/Ti at T = 8 K for different ap- 



FIG. 6. I/T2 as a function of temperature. Inset: Spin- 
echo decays are plotted as a function of r at two different 
temperatures. The solid line is the fit using Eq. (5). 

plied fields in order to check the effect of spin diffu¬ 
sion. It is known that in Heisenberg magnets, spin cor¬ 
relation function exhibits a diffusive effect due to the 
contribution of long wavelength or small wave vector 
(q ~ 0).^® In one-dimensional (ID) spin systems, spin dif¬ 
fusion shows a field dependence of 1/Ti oc l/\/ll type. 
Similar type of field dependence is reported in Heisen¬ 
berg spin-chain compounds (CH 3 ) 4 NMnCl 3 , AgVP 2 S 6 , 
CuCl2-2NC5H5, and SraCuOa.^e-ao Qn the other hand, 
in two-dimensional (2D) systems, spin diffusion results in 
a log(l///) dependence of 1/Ti, which has been exper¬ 
imentally observed in the quasi-2D Heisenberg magnet 
Zn 2 V 0 (P 04 ) 2 .^^’^^ As presented in the inset of Fig. 5(a), 
1/Ti for Li2CuW208 at T = 8 K is almost independent 
of the applied field confirming the absence of diffusive 
dynamics. 

The spin-spin relaxation rate I/T 2 was obtained by 
monitoring the decay of the transverse magnetization 
(Mxy) after a 7r/2 — t — tt pulse sequence as a function of 
the pulse separation time r and fitting by the following 
equation 

Mxy = -b C. (5) 

The extracted I/T 2 is plotted as a function of temper¬ 
ature in Fig. 6 and is found to be almost temperature 
independent. 

D. ^Li NMR spectra below Tn 

Field-sweep ^Li NMR spectra of Li2CuW208 below Tn 
are shown in Fig. 7. As discussed before, the NMR line 
above Tn broadens systematically with decreasing tem¬ 
perature. While approaching Tn, it broadens abruptly 
and then transforms into a nearly rectangular pattern 
well below Tn. The broadening of the spectra below Tn 
indicates that the Li site is experiencing the static in¬ 
ternal field in the ordered state through the hyperfine 
coupling between the ^Li nuclei and the ordered Cu^'*' 
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FIG. 7. Temperature-dependent field-sweep ^Li NMR spectra 
measured at 24.79 MHz below Tn. The red solid line is the 
fit to the spectra at 3.43 K using the sum of two Gaussian 
functions. 



FIG. 8. Temperature dependence of the internal field Hint 
obtained from ^Li NMR spectra measured at 24.79 MHz in 
the ordered state. The solid line is the fit by Eq. (6) as de¬ 
scribed in the text. Inset: Hnt vs. r and the solid line is the 
simulation of 0.053 x taking Tn — 3.51 K. 


moments. At low temperatures, the rectangular spectra 
develop two horns on either side of its edges. The spec¬ 
tral width is found to be increasing progressively upon 
lowering the temperature. 

The internal field Hint below Tn, which is proportional 
to the Cu^+ sublattice magnetization, was determined by 
taking half of the FWHM of the spectral lines. The tem¬ 
perature dependence of Hint after subtracting the high 
temperature (T > Tn) paramagnetic line width is shown 


in Fig. 8. In order to extract the critical exponent (/?) of 
the order parameter (sublattice magnetization), Hint(T) 
was fitted by the power law 

^fint(T) = Ho (^1 - . (6) 

It is noticed that Hint decreases sharply on approaching 
Tn. To get a reliable estimation of the critical exponent 
/3, one needs more data points close to Tn. Therefore, 
we have measured spectra in small temperature steps 
just below Tn (i.e. in the critical region). Our fit by 
Eq. (6) in the temperature range from 3 K to 3.6 K (see 
Fig. 8) yields /3 = 0.35(2) with Hq = 0.053(2) T and Tn 
= 3.51(1) K. In order to magnify the fit in the critical 
region. Hint is plotted against the reduced temperature 
T = in the inset of Fig. 8. The solid line is the fit 

by 0.053 X where Tn is fixed to the value 3.51 K. At 
low temperatures, the Hint value almost reaches satura¬ 
tion but much faster than the mean-field prediction (see 
the deviation between hollow spheres and solid line in 
Fig. 8). This behavior resembles temperature evolution 
of the ordered magnetic moment determined by neutron 
diffraction.^^ However, the neutron data are less accurate 
in this case because of difficulties in measuring very weak 
magnetic reflections in the vicinity of Tn, especially in a 
powder experiment. 

The value of /3 obtained from the analysis of sub-lattice 
magnetization reflects the nature of the critical fluctu¬ 
ations and hence the universality class or spin dimen¬ 
sionality of the spin systems. For a purely 3D ordering, 
the correlations in the critical regime, very close to Tn, 
will reflect 3D character. A tiny in-plane anisotropy can 
dominate the fluctuations if the dominant AF correlation 
length is large. If such a 3D ordering is still governed by 
the 2D processes, the phase boundaries that are observed 
should reflect the underlying 2D character. For instance, 
in BaNi 2 (P 04 ) 2 , P — 0.33 close to 3D model was found 
near Tn and another exponent /3 ~ 0.23 close to that of 
the 2D XY model was found in a broader region below 
Tn.^^ However, in some cases this 3D critical region is 
over such a narrow temperature range that it is not ac¬ 
cessible experimentally. Such a scenario has been realized 
before in the frustrated AFM square-lattice compounds 
Pb 2 V 0 (P 04)2 and Li 2 V 0 Si 04 . In these compounds, 
neutron diffraction experiments show that the ordering 
is 3D. In contrast, a reduced value of /3 close to the value 
of the 2D XY model has been reported from NMR mea¬ 
surements and the authors suggested that the transition 
to the columnar phase might be driven by the XY type 
anisotropy.^"^’^® The /? values expected for different uni¬ 
versality classes (spin and spatial dimensionalities) are 
tabulated in Ref. 34. Our obtained value of /3 ~ 0.35(2) is 
close to the value expected for any of the 3D spin models 
(Heisenberg, Ising or XY). This also reflects dominant 3D 
correlations in the critical regime and is consistent with 
our neutron diffraction experiments that predict the 3D 
ordering below Tn. 
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FIG. 9. (Color online) x{T) of Zn-doped Li2CuW208 samples 
measured at fioH = IT. Inset: Curie-Weiss temperature 6cw 
as a function of doping concentration x and the solid line is a 
linear fit. 


E. Zn doping 



FIG. 10. (Color online) Cp/T of Zn-doped Li2CuW208 sam¬ 
ples measured in zero magnetic field. The downward arrows 
point to the Neel temperatures Tjv. The solid lines are the 
fits by power law below Tn. 


As mentioned in Sec. II, all Zn-doped samples crys¬ 
tallize in triclinic symmetry, similar to the parent com¬ 
pound. The sample with x = 0.30 contained trace 
amounts of impurity phases, so its doping level may be 
slightly below 30%, but this minor deviation did not have 
any visible effect on the results. 

x(T) of doped samples normalized to one mole of Cu^'*' 
spins is shown in Fig. 9. The susceptibility maximum 
is systematically shifted to higher values of y and al¬ 
most temperature-independent, although it is gradually 
smeared out by the low-temperature susceptibility up¬ 
turn that becomes more pronounced at higher doping 
levels because of the increasing amount of paramagnetic 
spins. The weak kink at Tn vanishes upon doping. How¬ 
ever, the magnetic transition is still clearly seen in the 
heat-capaciy data presented below. 

The high-temperature data, x{T), was fitted by 


X{T) = Xo + 


C 

T + 0cw’ 


( 7 ) 


where xo is the temperature-independent susceptibility 
and the second term is the Curie-Weiss (CW) law with 
C being the CW constant and dew the CW temperature. 
From the fit above 50 K, the value of C was found to vary 
from 0.362 cm^ K/mol Cu to 0.454 cm^ K/mol Cu which 
correspond to the effective moment of /Teg = 1.7—1.9 /iu. 
The extracted values of dew decrease linearly upon dop¬ 
ing (see the inset of Fig. 9) from 12 K for the parent 
compound to 7.3 K for the x = 0.30 sample. This re¬ 
duction in dew confirms the dilution of the spin lat¬ 
tice. For random dilution, one expects that the slope 
of dcw(a^)/^cw(a^ = 0) vs. T is 1.0, which is comparable 
to the slope of 1.3 extracted from our data. 

For the sake of better presentation, we use a different 
scaling for the heat capacity and normalize the data to 
one mole of the compound. Figure 10 presents Cp/T vs. 
T for all the doped samples. The systematic reduction in 
the heat capacity maximum around 6 K follows the re¬ 


duced amount of the magnetic Cu^'*" ions. The position 
of the broad specific-heat maximum, which is only par¬ 
tially visible at a: = 0 and becomes very pronounced at 
X = 0.3, is roughly unchanged up to the highest possible 
doping level. 

Magnetic ordering transition is clearly seen in the heat- 
capacity data. The Neel temperature determined with 
the 0.05 K uncertainty from the maximum of the transi¬ 
tion anomaly, displays a systematic reduction from 3.9 K 
in the parent compound to 530 mK at a: = 0.25 and even¬ 
tually below 350 mK at a; = 0.30. Below Tn, Cp(T) data 
were fitted by the power law Cp oc T“ for all the doped 
samples, where a is the exponent. The a parameter is 
also strongly dependent on the doping level and decreases 
from Q; = 3ata: = 0toa~lata; = 0.25 (see Fig. 12). 


IV. DISCUSSION 
A. Static properties 

The broadening of the spectra below Tn suggests that 
a net field exists at the Li-site due to the nearest neigh¬ 
bor Cu2+ spins. As demonstrated in Fig. 1, Li atoms 
are located slightly off the plane (at an angle ~ 32°) but 
equidistant from two neighboring spins and each Li is 
coupled to two neighboring spins having the same direc¬ 
tion (either both up or both down). In such a scenario, 
one expects a strong internal field since the hyperfine 
fields due to the two up- or down- spins will be added 
up at the Li site. This is possible only if Li is coupled 
strongly to the Cu^’*' ions. As shown in Fig. 8, the in¬ 
ternal field at 1.8 K is about 0.038 T, which is relatively 
weak given the fact that Li is coupled to two spins hav¬ 
ing the same direction and, thus producing same direc¬ 
tion of the local field at the Li site. This weak internal 
field is comparable to the value reported in NMR for 














Zn 2 V 0 (P 04 ) 2 , where P is strongly coupled to one up- 
and one down- spins, so that the resulting internal field is 
only due to the difference in the hyperfine couplings along 
different pathways.Such a comparison proves 

that even though Li is coupled to two spins having the 
same direction, the diminutively small transferred hyper¬ 
fine coupling results in only a weak static field. The weak 
hyperfine coupling is also confirmed by the temperature- 
independent line shift at high temperatures (T > Tn)- 
Our spectral width analysis indeed suggests that the cou¬ 
pling between Li and Cu atoms is dipolar in nature, and 
the overlap of the s-orbital wave function of Li with the 
Sd wave function of Cu^’*' is almost negligible. 

For a powder sample in the magnetically ordered state, 
the direction of the internal field idint is distributed ran¬ 
domly in all directions. When the external field (H) is ap¬ 
plied, the effective field iLgff at a nuclear site is idint +H. 
In a ferromagnetic sample, iLeff is usually parallel to H 
and has a unique value, idint -I- id or jidint — id | in its 
magnitude, if idint is unique. This is because the ferro¬ 
magnetic moments are rotated easily to the direction of 
id in a relatively weak field. On the other hand, in case of 
an antiferromagnet, the moments hardly change their di¬ 
rections in the external field in usual NMR experiments. 
The direction of Hint has the random distribution even 
under external fields. Thus, has a distribution in 
its direction and its magnitude varies from | idint — id | to 
idint +id. In such a case one expects a typical rectangular 
powder spectra in the commensurate AFM ordered state 
which has been experimentally observed in compounds 
like CUV 2 O 6 , BaCo 2 V 208 , and BiMn 2 P 06 .^®^^® 

Our ^Li NMR spectra however exhibit two peaks on 
the top of a rectangular background. Such a spectrum 
with two horns has often been observed in NMR exper¬ 
iments on powder samples. For instance, powder NMR 
spectra of YMn 2 show a double horn pattern in the or¬ 
dered state where Mn moments are coupled antiferromag- 
netically and collinearly.^° Recently, the powder spec¬ 
tra of Zn 2 V 0 (P 04)2 in the commensurate AFM ordered 
state showed similar pattern with two horns on top of a 
rectangular background.According to Yamada et 
the shape of the spectra is also sensitive to instrumental 
features. When the axis of the sample coil of the NMR 
probe head, which is usually set perpendicular to the ex¬ 
ternal field dd, is not perpendicular to the distributed di¬ 
rections of ddeff’s then the intensity of NMR signal from 
the nuclei under investigation is reduced. In this case, 
the spectrum is curved around the middle having two 
peaks on either sides. Thus the double horn pattern in 
Fig. 7 is reminiscent of a commensurate AFM order and 
is consistent with the neutron diffraction results. 

As found from the neutron diffraction experiments,^^ 
the ordered magnetic moment below T/v is reduced to 
0.65 /is from its classical value of 1 fiB which has been 
attributed to strong quantum fluctuations. Alternatively, 
it can be interpreted as a coexistence of ordered and dis¬ 
ordered phases below T/v- In such a case, a narrow cen¬ 
tral line corresponding to disordered spins would have 



FIG. 11. ^Li 1/Ti for Li 2 CuW 20 g vs. the reduced tempera¬ 
ture (r). The solid line is a fit by the power law, l/Ti oc r *' 
with 7 ~ 0.1 and a fixed Tn — 3.6 K. 


been observed at the zero-shift position on top of the 
rectangular back ground. The ratio of area under the 
narrow line with the area under the broad back ground 
would give information about the fraction of the disor¬ 
dered phase below Tn. Our experimental spectra below 
Tn don’t show any narrow central peak suggesting that 
the ordered phase is purely collinear below Tn, without 
any signature of disorder. 


B. Dynamic properties 


Generally, one can express in terms of dynamic 
susceptibility x((f, w) as^°’^^ 


1 


K 


EI ^(9) 


|2 x"(t^o) 

Wo 


( 8 ) 


where the sum is over the wave vectors q within the 
first Brillouin zone, is the form factor of the hy¬ 
perfine interactions as a function of q, and is 

the imaginary part of the dynamic susceptibility at the 
nuclear Larmor frequency wq. For q = 0 and wq = 0, 
the real component of corresponds to the uni¬ 

form static susceptibility x- Therefore, one would expect 
l/yTiT to be temperature-independent in the param¬ 
agnetic regime. As shown in Fig. 5(b), l/yTiT is al¬ 
most temperature-independent above 30 K as expected. 
For T < 30 K, it increases slowly and shows a sharp in¬ 
crease below 10 K. This increase in 1/yTiT implies the 
dominance of Xjvf(T‘^o) over x(0, 0), which is due to the 
growth of spin fluctuations with q ^ Q, even at tempera¬ 
tures much higher than Tn. So, we can infer that strong 
antiferromagnetic spin fluctuations persist in a wide tem¬ 
perature range above Tn. This indicates the frustrated 
nature of the Li 2 CuW 208 in agreement with thermody¬ 
namic measurements.^^ 
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Since the dipolar coupling between Li and Cu atoms 
is known, it is possible to calculate the leading exchange 
coupling between Cu^“'" ions from the high-temperature 
value of 1/Ti. At high temperatures, 1/Ti is constant 
and can be expressed as:^^ 

/^ {'Jn9IJ.b)'^V^z'S{S + 1 ) 

\7iy x^oo SoJex \ Z' J 

where to ex = (| •/maxi fee/fi) \/2zS{S + l)/3 is the Heisen¬ 
berg exchange frequency, 2 : is the number of nearest- 
neighbor spins of each Cu^“'" ion, and z' is the num¬ 
ber of nearest-neighbor Cu^“'" spins for a given Li site. 
The value of ~ 7.83 x 10^^ cm“^ corresponds to 
1.1 kOe//rB.^^ Using the relevant parameters, ~ 
lAkOe/nB, IN = 103.962 x 10®rad z = 2, 

z' = 2, g = 2.17 [obtained from the x(T) analysis in 
Ref. 13], S = and the high-temperature (45 K) re¬ 
laxation rate of ( 7 ^ ) ~ 113.37s“^ for the Li site, 

the magnitude of the maximum exchange coupling con¬ 
stant is calculated to be Jmax/fcB — 13 K, which is in 
reasonable agreement with our leading exchange cou¬ 
pling (J/fee = 15 K) estimated from the band structure 
calculations.^® 

As discussed before, 1/Ti decreases slowly for T < 
20 K in contrast to the increase expected due to 
the growth of antiferromagnetic correlations. Simi¬ 
lar kind of T-dependence has been observed in sev¬ 
eral antiferromagnets®'*’^^^^"' and the decrease in 1/Ti 
above Tn is explained by the partial filtering of anti¬ 
ferromagnetic fluctuations from the neighboring spins. 
In those compounds, the probed nucleus is coupled to 
nearest neighbor spins having opposite directions. In 
contrast, the Li nuclei in Li2CuW208 are coupled to 
two nearest neighbor spins having the same direction. 
In this case, the filtering of antiferromagnetic fluctua¬ 
tions should not happen. Similar behaviour has been 
observed in frustrated breathing pyrochlore lattice com¬ 
pound LiInCr408, where the decrease in 1/Ti is ex¬ 
plained by the opening of a spin gap above Tn.^® In 
LiInCr 408 , x('T) also shows an activated behaviour at 
low temperatures.^® However, in this compound xi'^) 
does not show any exponential decrease, thus ruling out 
the possibility of a spin gap at low temperatures. For 
spin-1/2 Heisenberg AFM spin chains, it has been pre¬ 
dicted theoretically that the contribution of the stag¬ 
gered component {q = zLir/a) to the spin-lattice relax¬ 
ation rate behaves as 1/Ti ~ T® while the contribution 
of the uniform component {q = 0) scales as 1/Ti ~ 

The contributions of staggered and uniform components 
normally dominate in the low-temperature (T ^ J) and 
high-temperature (T ~ J) regimes, respectively.^® Our 
experimentally observed linear behavior of 1/Ti with 
temperature over the 5.5 K < T < 16 K range in 
Fig. 5(a) is possibly due to ID physics. Indeed, we 
found moderate ID anisotropy of the 3D spin lattice in 
Li 2 CuW 208 with the leading interaction running along 
the 6-direction.®® However, the purely ID model is by far 


insufficient to describe the physics of the compound even 
in the T ~ J temperature range.®® 

At the 3D ordering temperature, the correlation length 
is expected to diverge, and 1 /Ti in a narrow temperature 
range just above Tn (i.e., in the critical regime) should 
be described by the power law, l/Ti oc t“^, where 7 is 
the critical exponent. The value of 7 represents universal 
property of the spin system depending upon its dimen¬ 
sionality, the symmetry of the spin lattice, and the type 
of interactions. To analyze this critical behavior, I/Ti 
is plotted against the reduced temperature r in Fig. 11. 
The data just above Tn (t < 0.5) were fitted by the 
power law with a fixed Tn — 3.6 K yielding 7 ~ 0.10. 
For the 3D Heisenberg antiferromagnet, a mean-field the¬ 
ory predicts 7 = ^ and a dynamic scaling theory gives 
7 = |.®®d® Our experimental value of 7 ~ 0.10 is far 
below these theoretically predicted values. This effect re¬ 
quires further investigation and may be intertwined with 
the strong frustration in 3D that affects spin dynamics 
above Tpf. 

In the magnetically ordered state (T < Tn)) the strong 
temperature dependence of 1/Ti is a clear signature of 
the relaxation due to scattering of magnons by the nu¬ 
clear spins.®® For T A/fee, l/7i follows either a T® 
behavior or a T® behavior due to a two-magnon Ra¬ 
man process or a three-magnon process, respectively, 
where A is the energy gap in the acoustic magnon 
spectrum.®®’®® For instance, in Pb 2 V 0 (P 04 ) 2 , 1/Ti fol¬ 
lows a T® behaviour®"® while a T® behaviour has been re¬ 
ported for spin-/ square-lattice compound Zn 2 V 0 (P 04)2 
(Ref. 32) and the decorated Shastry-Sutherland lattice 
in the spin-/ compound CdCu 2 (B 03)2 (Ref. 49). On 
the other hand, for T ^ A/Zcb; h follows an activated 
behavior 1/Ti (x 

Clearly the T® behaviour couldn’t describe our low 
temperature l/Ti below Tn for Li 2 CuW 208 . As one can 
see in Fig. 5(a), we fitted 1/Ti below Tn by a T® be¬ 
havior. An activated behavior would be consistent with 
our data too. However, the fit by 1/Ti oc 7 ’ 2 g-A/fcBT 
the temperature range 2 K to 3.3 K yields A ~ 3.7 K 
that is unusually large for this system and amounts to 
about 25% of Jmax- This substantial spin gap is hard 
to justify for our system, because inversion symmetry 
of the crystal structure eliminates any Dzyaloshinsky- 
Moriya couplings, which are the leading anisotropy term 
in Cu2+ oxides.®® The remaining anisotropy is of sym¬ 
metric (Ising) type and does not exceed 1 — 2 % in 
cuprates,®® so it can not be responsible for the spin gap 
extracted from the fit of 1/Ti. It is also worth noting that 
the specific heat follows the T® behavior at least down 
to 1.8 K (Fig. 10). Therefore, we conclude that the T® 
behavior related to three-magnon processes is the most 
plausible scenario of the low-temperature spin-lattice re¬ 
laxation, although measurements below 2 K could be use¬ 
ful to check whether the activated behavior might set in 
at very low temperatures corresponding to smaller values 
of the spin gap. 



10 



FIG. 12. Log Tn and a vs. x obtained from the C'p(T) analysis 
in Fig. 10 in the left and right y-axes, respectively. The solid 
lines are the linear fits. Inset: Tt^/J' vs. LJ' where the 
average chain length L — 1/x — l and Jid = 15 K is compared 
with the theoretical results from Ref. 54. 


C. Effect of doping 

The spin lattice of Li2CuW208 can be diluted by Zn 
atoms, and relatively high doping levels of 25 — 30 % can 
be reached, which is not very common for oxides 

because Zn^+ does not favor the four-fold square-planar 
oxygen environment, which is triggered by the Jahn- 
Teller effect for Cu^'*'. Other dopants may be compatible 
with this local environment of the Cu^+ site, but their 
size difference prevents high doping levels. For exam¬ 
ple, only 3 % of Pd^+ can be doped into the spin-chain 
antiferromagnet Sr 2 Cu 03 .®® 

The linear evolution of Oqw upon doping (see the in¬ 
set to Fig. 9) confirms that Zn^+ ions randomly break 
exchange bonds and reduce the overall exchange energy. 
This should reduce the ordering temperature Tn, whereas 
the slope of Tn vs. x depends on the nature of the sys¬ 
tem. Quantum systems, especially those having reduced 
dimensionality, will show higher slopes than their classi¬ 
cal counterparts because the effect of dilution is ampli¬ 
fied by quantum fluctuations that impede the long-range 
magnetic ordering. In spin-^ square-lattice antiferromag- 
nets, the slope of TN(a;)/TN(0) vs. x is expected to be 
around 3.2.®® Even higher values were reported exper¬ 
imentally and ascribed to the combined effect of dilu¬ 
tion, low-dimensionality and frustration induced by the 
doping.®^ Remarkably, none of those scenarios can be ap¬ 
plied to Li 2 CuW 208 , where the doping dependence of Tn 
can be tentatively linearized in the log Tn vs. x (and not 
in Tn vs. x) coordinates. This exponential reduction 
in Tn is strongly reminiscent of chain antiferromagnets. 
However, the model reported in the literature®'^ provides 
only a very crude fit of our data assuming a weak inter¬ 
chain coupling J'! JvD = 0.065 (see the inset of Fig. 12). 
Comparing our data with the theory reported in Ref. 58, 
we conclude that Li2CuW208 is similar to quasi-lD an¬ 
tiferromagnets with very weak interchain couplings, be¬ 


cause the crossover between the exponential and linear 
decay of Tn is seen at J' j Jro — 0.2. 

Naively, one would relate the ID-like evolution of Tn 
upon doping to the presence of broad maxima in the 
magnetic susceptibility and specific heat. However, the 
weakly curved magnetization isotherm®® of Li 2 CuW 208 
contradicts this ID scenario and indicates strong inter¬ 
chain couplings that, according to our ab initio results,®® 
reach 40 — 50 % of Jid rendering Li2CuW208 a spa¬ 
tially anisotropic 3D antiferromagnet rather than a spin- 
chain antiferromagnet. Moreover, the ordered moment 
of 0.65 fiB (Ref. 13) also indicates®® interchain couplings 
of J'/Jm = 0.4 — 0.5 that are much stronger than 
J'/Jm = 0.065 expected from the doping dependence 
ofTN. 

Altogether, the doping experiments reveal peculiarities 
of Li2CuW208. The evolution of Tn in doped samples 
gives insight into the nature of the parent compound, 
where strong magnetic frustration in 3D affects the or¬ 
dered magnetic moment in the ground state®® and fa¬ 
cilitates the suppression of Tn upon the dilution. The 
ground state of doped samples may be quite interesting 
on its own, and we note in passing that the power-law 
dependence of the low-temperature specific heat changes 
systematically from the conventional T® behavior in the 
parent compound to the nearly linear behavior at a; = 
0.25. This effect warrants further investigation that, 
however, goes beyond the scope of our present study. 


V. SUMMARY 

Frustrated nature of the 3D antiferromagnet 
Li2CuW208 has been probed by ''Li NMR mea¬ 
surements and experiments on magnetically diluted 
Zn-doped samples. The double-horn spectral shape 
of the NMR line below Tn is attributed to the for¬ 
mation of commensurate magnetic order and the 
non-perpendicularity between the axis of the NMR coil 
and the resonance field. The absence of a narrow central 
line at the zero-shift position also supports earlier 
neutron diffraction results suggesting commensurate 
magnetic order with the reduced ordered moment. In 
contrast to the isostructural compounds Li 2 NiW 208 and 
Li 2 CoW 208 , no magnetic transitions beyond that at Tn 
have been observed. 

The nature of the 1 /Ti process below Tn remains some¬ 
what ambiguous and a better analysis requires more data 
points in a wider temperature range below 2 K. The anal¬ 
ysis of 1/xTiT suggests the frustrated nature of the com¬ 
pound. The critical exponent (3 obtained from the sub¬ 
lattice magnetization below Tn is consistent with the 3D 
nature of the spin lattice. 

In Zn-doped sample, Tn decreases exponentially upon 
non-magnetic dilution x. This behavior is reminiscent 
of quasi-ID antiferromagnets and may be a footprint of 
magnetic frustration in the parent compound, because 
the quasi-lD scenario is clearly excluded by the size of the 
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ordered moment and by the microscopic analysis. Heat 
capacity below Tn for all the doped samples follows the 
power law (Cp oc T“), but the exponent (a) decreases 
from 3 to 1 as a: increases from 0 to 0.25. 
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